Evolution of the cosmic ray anisotropy above 10 14 eV 
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ABSTRACT 

The amplitude and phase of the cosmic ray anisotropy are well established experimentally 
between 10 11 eV and 10 14 eV. The study of their evolution into the energy region 10 14 — 10 16 
eV can provide a significant tool for the understanding of the steepening ( "knee" ) of the primary 
spectrum. In this letter we extend the EAS-TOP measurement performed at Eq rs 10 14 eV, to 
higher energies by using the full data set (8 years of data taking). Results derived at about 10 14 
and 4 • 10 14 eV are compared and discussed. Hints of increasing amplitude and change of phase 
above 10 14 eV are reported. The significance of the observation for the understanding of cosmic 
ray propagation is discussed. 



Subject headings: Cosmic rays, diffusion 

1. Introduction 

The steepening ("knee") observed at E ~ 3 • 
10 15 eV represents a main feature of the energy 
spectrum of cosmic rays and its characterization 
is therefore a main tool for the understanding of 
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the galactic radiation. Composition studies have 
shown that it is related to the steepening of the 
lightest primaries (protons, helium, CNO) spectra 



(jAglietta et al. Il2004 lAntoni et al. 1120051 ) 



Such effect can be due, on the one side, to 
energy limits of the acceleration process at the 
source, namely diffusive shock acceleration in 
supernova remnants, generally considered to be 
the sources of galactic cosmic rays. The max- 
imum energy of the accelerated protons is, in- 
deed, c alculated to occur in the 10 15 eV energy 
regio n 
l2007h. but 



([Berezhko et al I Il996t berezhko fc Volk 



could reach up 



to about 10 17 eV 
On the other side, 



(jPtuskin fc Zirakashvili 1120031) 
this feature has been possibly explained in terms of 
a change in the cos mic ray propa g ation properties 
insid e the Galaxy (jPeters 1960; IZatsepin et al. I 
1962). Galactic propagation is described through 



diffusion models whose parameters have been ob- 
tained through composition studies (mainly from 
the ratio of secondary to primary nuclei) at en - 
ergies well below 1 TeV (see e.g. (jJones et al 
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20011 IStrong et al. I [20071 ) ). The diffusion coeffi- 
cient, D, is found to increase with magnetic rigid- 
ity (D oc R 6 , or D oc i? 3 for models includ- 
ing reacceleration) . However, no confirmation, 
and no information has till now been obtained 
at higher energies, where the main observable is 
represented by the large scale anisotropy in the 
cosmic rays arrival directions, that is known to 
be st rictly related to th e diffusion coefficient (see 



e.g. ( Berezinskv et al. 1990l) ). The study of the 



evolution of the anisotropy in the "knee" energy 
region can therefore provide a significant test of 
the diffusion models, and a valuable insight for the 
discrimination between the two possible explana- 
tions of the spectral steepening. 

At En « 10 14 eV the EAS-TOF0 results 
(jAglietta et al. I Il996h demonstrated that the 
main features of the anisotropy (i.e. of cosmic 
ray propagation) are similar to the ones measured 
at lower energies (10 11 10 14 eV), both with re- 
spect to amplitude ((3 -j- 6) • 10~ 4 ) and phase ((0 



4) h LST) dGombosi et al. lll975UFenton et al. 


1975 


: iNaeashima et al. 


1989; Alekseenko et al. 


1981 


: Andreev et al. Hl987l:lAmbrosio et al. 112003: 


Munakata et al. 1997: 


Amenomori et al. 12005; 


Guillian et al. 2007; Abdo et al. 20081). At 



higher energies the limited statistic s does not al- 



2007) 



low to draw any firm conc l usion (IKifune et al 
1984t iGherardv et al. I [1981 lAntoni et al. 1 120041 



Amenomori et al. 112006c IQver et al. 

In this letter we present the EAS-TOP mea- 
surement based on the full data-set and we extend 
the analysis to about 4 • 10 14 eV. 

2. The experiment and the analysis 

The EAS-TOP Extensive Air Shower array was 
located at Campo Imperatore (2005 m a.s.L, lat. 
42° 27' N, long. 13° 34' E, INFN Gran Sasso Na- 
tional Laboratory). The electro magnetic detec 
tor (u sed for the present analysis) (jAglietta et al 



1993) consisted of 35 modules of scintillator coun- 
ters, 10 m 2 each, distributed over an area of about 
10 5 m 2 . The trigger was provided by the coinci- 
dence of any four neighbouring modules (threshold 
n p k 0.3 m.i. p. /module), the event rate being / « 
25 Hz. The data under discussion have been col- 
lected between January 1992 and December 1999 



1 The Extensive Air Shower array on TOP of the Gran Sasso 
underground laboratories. 



for a total of 1431 full days of operation. 

To select different primary energies, a cut 
is applied to the events based on the number 
of triggered modules (see table 1). The av- 
erage primary energies are evaluated for pri- 
mary proto ns and QGSJET01 h adron interac- 



tion model (IKalmykov et al. 1119971) in CORSIKA 



dHeck fc Knapp 1119981) . 



Table 1: Characteristics of the two classes of 
events used in the analysis: number of triggered 
modules, primary energy and number of collected 
events in the East+ West sectors. 



Class 


N modules 


E 


[eV] 


New 


I 


> 4 


1.1 


10 i4 


1.5 • 10 y 


II 


> 12 


3.7 


10 14 


1.7- 10 s 



For the analysis of the anisotropy, we adopt a 
method based on the counting rate differences be- 
tween East-ward and West-ward directions, that 
allows to remove counting rate variations of at- 
mospheric origin. The events used in the anal- 
ysis (see table 1) are the ones with azimuth an- 
gle <j> inside ±45° around the East and West di- 
rections, and zenith angle 8 < 40°. The differ- 
ence between the number of counts measured from 
the East sector, Ce{P), and from the West one, 
Cw(t), at time t in a fixed interval (At = 20 
min), is related to the first derivative of the in- 
tensity I(t) as: § ~ D(t) = g gM-CV(t) where 



st 



5t is the average hour angle between the vertical 
and each of the two sectors (1.7 h in our case). 
The harmonic analysis is performed on the differ- 
ences D(t); the amplitudes and phases of the vari- 
ation of I(t) are obtained through the integration 
of the corresponding terms of the Fourier series 



(jAglietta et al. 1120071 ) 



3. Results 

The harmonic analysis has been performed in 
solar, sidereal and anti-sidereal timed. We de- 
scribe in sect 3.1 the results of the analysis, while 



2 The anti-sidereal time is a fictitious time scale symmetrical 
to the sidereal one with res pect to the solar time a nd that 
reflects seasonal influences (Farley & Storey 1954). 
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Eo{eV] 


A L io 4 <t>l iM pU%) 


A L io 4 tijh] pU%) 


< si ,i io 4 <f>LM pLM 


1.1 ■ 10 14 
3.7- 10 14 


2.8 ±0.8 6.0 ±1.1 0.2 
3.2 ±2.5 6.0 ±3.4 44.1 


2.6 ±0.8 0.4 ±1.2 0.5 
6.4 ±2.5 13.6 ±1.5 3.8 


1.2 ±0.8 23.9 ± 2.8 32.5 
3.4 ± 2.5 22.3 ± 3.2 39.7 




Mil io 4 &] pL t i(%) 


a& io 4 <j>llM P&{%) 


Km io 4 KiidM pIU%) 


1.1 • 10 14 
3.7- 10 14 


1.4 ±0.8 7.0 ±1.2 21.6 
1.7 ±2.5 - 79.4 


2.3 ±0.8 6.3 ±0.7 1.6 
1.5 ±2.5 - 83.5 


0.6 ±0.8 - 75.5 
1.2 ±2.5 - 89.1 



Table 2: Results of the analysis of the first (amplitude A 1 , phase (f) 1 , and Rayleigh imitation probability P 1 ) 
and second harmonic (A 11 , (f) 11 , P 11 ) in solar (columns 2-4), sidereal (columns 5-7), and anti-sidereal time 
(columns 8-10). Phases are not defined when amplitudes are smaller than their uncertainties. 



in sect. 3.2 we show the related counting rate 
curves. 

3.1. The harmonic analysis 

For the two different primary energies, the re- 
constructed amplitudes and phases of the first and 
second harmonics are shown in table 2, together 
with the corresponding Rayleigh imitation proba- 
bilities (P). 

- Concerning the first harmonic: 

(a) At 1.1 • 10 eV, from the analysis in 
solar time , the obtained amplitude and phase 
{A[ ol = (2.8 ± 0.8) ■ IO" 4 , 4>L = (6.0 ± 1.1) h, 
Psoi = 0.2%) are in excellent agreement with the 
expected ones from the Co mpton-Getting effect 



(jCompton fc Getting I 119351 ) due to the revolu- 
tion of the Earth around the Sun: at our latitude 
^ oi , CG =3.0-10- 4 , so ;,cg=6.O h. 

With respect to the sidereal time analysis, the 
measured amplitude and phase {Al id = (2.6 ± 
0.8) ■ 10~ 4 , 4>i ld = (0.4 ± 1.2) h LST), with imita- 
tion probability = 0.5%, confirm th e previous 
EAS-TOP result jAglietta et al~lll996h . 



Bi-monthly vectors representing the first har- 
monic are shown in figure [T] (dots), together with 
the expected ones (stars) from the measured so- 
lar and sidereal amplitudes. The expected anti- 
clockwise rotation of the vector is clearly visible, 
showing that, at any time, the composition of the 
two vectors is observed, and that the expected and 
measured individual values are fully compatible 
within the statistical uncertainties. 

(b) At 3.7 ■ 10 14 eV the amplitude and phase 
of the measured first harmonic in solar time are 



still consistent with the expected ones for the so- 
lar Compton-Getting effect, although, due to the 
reduced statistics, the chance imitation probabil- 
ity is rather high. 




12 h 



Fig. 1. — Bi-monthly solar vectors at 1.1 • 10 14 eV. 
Capital and small letters refer respectively to the 
expected and theoretical points (A, a = Jan + Feb; 
B,b = Mar + Apr; C,c — May + Jun; D,d = Jul 
+ Aug; E,e = Sep + Oct; F,f = Nov + Dec). The 
typical vector statistical uncertainty is shown. 



Concerning the analysis in sidereal time , we ob- 



tain A{ ld = (6.4±2.5)-10" 



(13.6±1.5) h 



LST, with an imitation probability of about 3.8%. 
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This indicates therefore a change of phase (from 
0.4 to 13.6 h) and an increase of amplitude (by a 
factor 2.5) with respect to the first harmonic mea- 
sured at 1.1 • 10 14 eV. 

- Concerning the second harmonic most sig- 
nificant (P*L — 1.6%) is the amplitude observed 
in sidereal time in the lower energy class of events 
(comparable with the first harmonic one: = 
(2.3 ± 0.8) ■ 1CT 4 , d>{L = (6.3 ± 0.7) h LST) (see 



also (Alekse enko et al. II1981D ). 



Both at 1.1 ■ 10 14 eV and 3.7- 10 14 eV, no signif- 
icant amplitude is observed in anti-sidereal time , 
showing that no additional correction is required 
due to residual seasonal effects. 



3.2. The counting rate curves 
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Fig. 2. — Thick black lines: counting rate curves in 
solar (a), sidereal (b) and anti-sidereal (c) time at 
1.1 • 10 14 TeV. The statistical uncertainty for each 
bin is given in the first one. The curves resulting 
from the first harmonic analysis are also shown 
(light black lines); for the sidereal time curve, the 
combination of first and second harmonic ( dotted 
black line) is additionally superimposed. 



Besides the harmonic analysis, it is interesting 
to visualize the variations of the cosmic ray inten- 
sity versus time, I(t), as reconstructed by integra- 
tion of the East- West differences, D{t). They are 
shown in figs. [2] and [31 for the classes of events at 
1.1 • 10 14 eV and 3.7 • 10 14 eV, respectively (a,b,c 
for solar, sidereal, and anti-sidereal time scales). 

As already shown by the harmonic analysis, at 
both energies the curves in solar time are domi- 
nated by the Compton-Getting effect due to the 
motion of the Earth, and no modulation is visible 
in the anti-sidereal time scale. 
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A main difference is observed in the sidereal 
time curves: while the shape of the curve at 1.1 ■ 
10 14 eV is in remarkable agreement with the EAS 
and muon measurements reported at and below 
100 TeV, the curve related to the highest energy 
class of events is characterized by a broad excess 
around 13-16 h LST. 
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Fig. 3. — Thick black lines: counting rate curves 
in solar (a), sidereal (b) and anti-sidereal (c) time 
at 3.7 • 10 14 TeV. The curves resulting from the 
first harmonic analysis are also shown (light black 
lines). 



4. Conclusions 

High stability data obtained from long time 
observations (8 years) from the EAS-TOP ar- 
ray confirm the amplitude and phase of the cos- 
mic ray anisotropy already reported at 10 14 eV: 
Kid = (2-6 ±0.8) ■ 10- 4 , 4, = (0.4 ± 1.2) h LST, 
with Rayleigh imitation probability P/- d = 0.5%. 
The result is supported by the observation of the 
Compton-Getting effect due to the revolution of 
the Earth around the Sun, and by the absence of 



anti-sidereal effects. It confirms the homogeneity 
of the anisotropy data over the energy range 10 11 - 
10 14 eV. 

At higher energies (around 4 • 10 14 eV) the 
observed anisotropy shows a larger amplitude, 
Al ld = (6.4 ± 2.5) • 10~ 4 , and a different phase, 
<\>\ id = (13.6 ± 1.5) h LST, with an imitation prob- 
ability of 3.8%. The statistical significance is still 
limited, but the measurement has the highest sen- 
sitivity with respect to previous experiments at 
these energies, and it is not in contradiction with 
any of them. 

The dependence of the anisotropy amplitude 
over primary energy (A cx Eq) deduced from the 
present two measurements can be represented by 
a value of r5 = 0.74 ±0.41. Therefore, at least in 
the energy range (1 — 4) • 10 14 eV, such dependence 
is compatible with that of the diffusion coefficient 
as derived by composition measurements at lower 
energies. 

On another side, the sharp increase of the 
anisotropy above 10 14 eV may be indicative of a 
sharp evolution of the propagation properties, and 
therefore of the diffusion coefficient just approach- 
ing the steepening of the primary spectrum. This 
opens the problems of obtaining an improved the- 
oretical and experimental description of the whole 
evolution of the diffusion processes vs primary en- 
ergy, and understanding how such evolution could 
affect the energy spectra at the "knee". From the 
experimental point of view, the extension of the 
anisotropy measurements with high sensitivity to 
and above 10 15 eV will be of crucial significance. 
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